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In this paper, we studied a generalized Bose-Hubbard model on a checkerboard lattice with 
topologically nontrivial flat-band. We used mean- field method to decouple the model Hamiltonian 
and obtained phase diagram by Landau theory of second-order phase transition. We further calculate 
the energy gap and the dispersion of quasi-particle or quasi-hole in Mott insulator state and found 
that in strong interaction limit the quasi-particles or the quasi-holes also have flat bands. 



I. INTRODUCTION 

The Mott insulator (MI) - superfluid transition (SF) of 
ultracold bosons in optical lattices has become a hot 
topic in quantum simulation |3| and a great deal of works 
have been done. The Bose-Hubbard model proposed in 
Ref. [4] have been widely investigated both theoretically 
and experimentally so that a comprehensive understand- 
ing of this many-body system can be achieved. 

On the other hand, the states with topological or- 
der become another hot topic and have been explored 
from different aspects. Since first discovered in two- 
dimensional electron gas with Landau levels in strong 
magnetic field, integer and fractional quantum Hall ef- 
fect (IQHE and FQHE) have attracted great attention. 
Haldane pointed out that IQHE may appear in honey- 
comb lattice without Landau levels [5|- However, FQHE 
in a lattice model without Landau levels had not been 
discovered until the proposal of the lattice models with 
topologically nontrivial flat-band based on the mecha- 
nism of quadratic band touching 6 8]. In these lattice 
models (kagome lattice, 7r-flux square lattice and hon- 
eycomb lattice), the hoppings of next-nearest or next- 
next-nearest (NNN and NNNN) neighbor besides nearest 
neighbor (NN) are introduced to achieve the flat-bands 
with a high flatness ratio (the flatness ratio defined as 
the band gap over the band width) in different lattice 
models. To obtain FQHE states, people always consider 
these tight-binding models with the NN and NNN repul- 
sions for fermions and corresponding tight-binding mod- 
els under a hard-core condition for bosons. Then the 
numerical results confirm the existence of topologically 
nontrivial flat-band. For example, in Ref. ([9]), in an ex- 
tended bosonic Haldane model with the NN and NNN 
repulsions under a hard-core condition, the v = 1/2 and 
the v = 1/4 bosonic FQHE states are predicted. 

So we may ask questions : in a bosonic lattice model of 
topologically nontrivial flat-band with on-site Coulomb 
interaction, what's the ground state and whether the 
topologically nontrivial flat-band affects the superfluid 
state and Mott state? This becomes the starting point 
of this paper. In this paper, we study the MI-SF transi- 



tion in a bosonic system of the checkerboard model with 
topological flat-band of which there is no NN or NNN 
repulsions. Instead, we replace the hard-core condition 
by a tunable on-site repulsive interaction and consider a 
generalized Bose-Hubbard model with the hopping pa- 
rameters as that of the models with topologically non- 
trivial fiat-band. We will use the mean field approach to 
obtain the quantum phase transitions between MI phase 
and SF phase [lCj, LLJJ. Then the elementary excitations 
are studied in Mott insulator state. We derive the energy 
gap and the dispersion of quasi-particle or quasi-hole. 

The paper is organized as follow: We first introduce 
the generalized Bose-Hubbard model with the hopping 
parameters as that of the models with topologically non- 
trivial flat-band in Sec. II, then study this model in Sec. 
Ill by mean field approach. And in Sec. Ill we obtain 
the phase diagrams of this model. The properties of col- 
lective modes are analyzed in the last subsection in Sec. 
IV. Finally we draw the conclusion in Sec. V. 



II. THE MODEL HAMILTONIAN 




FIG. 1: The scheme of the square lattice, of which the two 
different lattice sites in blue or green colors are two sublattices 
A and B. The NN hoppings t are represented by lines with 
arrows whose directions are sign of the phase ip. The dot- 
ted lines represent the NNN hoppings which are — t' between 
green sites (sublatt ice A) and t' between blue sites (sublatt ice 
B). Two dashed lines representing the NNNN hoppings are 
displayed in the scheme, too. 
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We first generalize the fermionic topological flat-band 
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model to bosonic system with the Hamiltonian 
H =~ t Yl [^^exp^)] ±t f h is 

(hi) {{3,1)) 

«<i,0» J 
-/i^6t.6 r . +H.c. 

3 

where b\ (b rj ) is the bosonic creation (annihilation) op- 
erator of the j-th site with position vector rj. 

and denote that j and I are the NN, 

NNN and NNN sites, respectively. The phase factor of 
NN hopping ipji equals ±(p (Figured)), £, t' ', t" are the 
hopping parameters for NN hoppings, NNN hoppings, 
NNNN hoppings, respectively. U is the on-site interac- 
tion strength and \i is the chemical potential. 

In the following parts we will use mean field approach 
to calculate above bosonic topological flat-band model. 
In addition, to show the effect of the flat-band on MI- 
SF transition, we also calculate the generalized Bose- 
Hubbard model without topological flat-band and com- 
pare their phase diagram. 



III. MEAN-FIELD APPROXIMATION AND 
PHASE DIAGRAM 

In this section we derive the decoupled effective Hamil- 
tonian from Eq. (pQ) by mean-field approximation and cal- 
culate the corresponding perturbation energy up to sec- 
ond order. Then the equation of critical line of MI-SF 
transition is obtained by Landau theory of phase transi- 
tion. 



A. Mean-field approximation 

In the strong-coupling Limit (t <C U), a localized su- 
perfluid order parameter ifj rj = (bj ) = (b rj ) is intro- 
duced so as to rewrite the hopping terms by a consistent 
mean-field method as 

h l 3 K - + b ri ip rj - ij) r . ijj ri . (2) 

Thus the Hamiltonian can be written as 



IT* =-tcos^Y,H 



3 1=1 

z' 



3 1=1 



z 



u 



- n rj (n rj - 1) - n ti, b rj , 



(3) 



where z, z f , z" denote the coordination numbers of NN, 
NNN and NNNN hoppings and r h rj, t" are the cor- 
responding translation vectors, respectively. Hence the 
decoupled effective Hamiltonian can be rewritten with 
respect to Yj as 

H*f = -zt cos <p(bl d + b r . - ipTjWrj+T 

-zt"(bi 3 +b Y3 -^ 3 )^ Y3+T „ 



u < 

+ y^K - V)-[in T , . 



(4) 



Since there are two sublattices of this lattice (A and 
B), we can define two order parameters ^a = (H jeA ) ano ^ 
ipB = (bt jeB ) that correspond to the Bose condensations 
on A site and on B site, respectively. The order param- 
eters satisfy the condition i^a = = ^yj+t' — ^Yj+t" 
and i(jb = ^Yj+t if site j denotes an A site and similar 
condition for B site can be obtained if we substitute A 
with B. Then we can construct an effective Hamiltonian 
of a two-site cell with the condition z = z' = z" as 

JT ff = (ztyHH! + H 2 ) = (zt)-\H*f + ff r f +T ) 

= 7^ ( n l + n 2 — n ) ~ M n + 2^1 ^2 COS if 



+ t"^i + t"^ A 2 + (^2 cosy? + *>i)Vi 

+ (^iCOS^ + ^Vs)^ , 



(5) 



where ^1 = ip A , ^2 = and U = ^, p, = f v t" = £, 
Vi = — (bj + bi). We then divide the effective Hamiltonian 
into two terms: the unperturbed term 



# (0) = ^ (n\ + nl - n)- fin + 2ip!ip 2 cos <p + t'Vj+t'V| , 

(6) 

and the perturbation term 

H' = (^i" + ^2 cos <p)Vi + (^1 cos <p + ^2t")V 2 . (7) 
Hence the unperturbed ground energy is given by 

£ (0) = ^ (n\ + n 2 2 -n)- fin + 2^2 cos p + FVi+^Vi • 

(8) 

As first-order perturbation of energy vanishes, we get 
second-order energy perturbation 



{j) 2 COS if + t"jj 1 ) 2 (U + 

'{-mU + p) [(m-VjU-ji] 
(V>i cos <p + i"ip 2 ) 2 (U + p) 



(9) 



(-n 2 U + fi) [(n 2 - 1)U - p] ^ 
Thus energy up to second-order is obtained as 

£(^1,^2) « 00 + 02^ + 02^1^2 +b 2 ^ + (D(i> 2 ) , (10) 
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where 



U 



cio =2"(»ii + n 2 — n) — fin 



a 2 =t" 



i" 2 (U + p) 



b 2 =i" 



c 2 =2 cos ip 



(-mU + p) [(ni - 1)U - p] 
(U + p) cos 2 ip 
(-n 2 U + p) [(n 2 -l)U-p] 

(U + p) cos 2 ip 
~ (-nJJ + p,) [(m - l)U-p] 

i" 2 (U + p) 
(-n 2 U + p) [(n 2 - 1)U - p] 
2i"(U + p) cosip 



(-mU + p) [( ni -i)U-p\ 
2i"(U + p) cosip 



(-n 2 U + p) [(n 2 - 1)U - p] 



(11) 



B. Phase diagram 

The MI-SF phase transition occurs on the condition 
when Gaussian curvature is zero at the point ipi — 
0, ^2 = 0, namely 
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It implies that 



4a 2 b 2 - <? 2 = . 



(12) 



(13) 



With m = n 2 = n, considering the symmetry between 
two sublattices, a 2 = b 2l we have 2a 2 = ±C2. Solve the 
equation for /i, we can get the boundary condition 



1 

2 

± 



(2n- l)/7- (£" - cosy?) 



Z7 2 



2(2n + - cos ip)U + (t" - cosip) 2 



(14) 



for negative U , namely positive U as t = — 1. 

The phase diagram is displayed in Figure [2j In Figure 
[2j there are three MI-SF phase transition lines for each 
particle number configuration. The particle number con- 
figuration {ni 577-2} represents that there are ri\ bosons 
on each A lattice site while n 2 bosons on each B lattice 
site. Due to the existence of inversion symmetry, we have 
n\ n 2 . 

The solid lines are phase boundaries of the flat-band 
case with t = — 1, i! — — - 1 7= , i" = —^-7= and ip = 

2+^2' 2+2^2 ^ 4_ 

We can see that quantum phase transition occurs when U 
increases by fixing ft. The phase diagram is similar to the 




FIG. 2: The phase diagram of MI-SF transition. The solid 
lines are the phase boundaries of topological non-trivial flat- 
band case with t — —1, F — — ? i" = 2+ 2^/2 an< ^ ^ ~ f • 
The dashed lines are the phase boundaries of the case with 
t = -1, ? = 0, t" = and y? = f . The dotted lines are 
phase boundaries of the case with t = — 1, t' - 
1 _ 

4 • 



and 



" 2(2+^2) ' 

On the other hand, the red lines 



2(2+2^) 

represent phase boundaries of configuration {1; 1}, blue ones 
{2; 2}, and green ones {3; 3}. 



traditional Bose-Hubbard model with only nearest hop- 
ping term[T3,[ll|. To show the effect of flat-band, we also 
calculate other two cases without topological flat-bands: 
the (^-flux Bose-Hubbard model with only nearest hop- 
ping term or i! = t" = 0, ip = j (the dashed lines) and 
the generalized Bose-Hubbard model with smaller NNN 
and NNNN hopping parameters compared with flat-band 



case, namely ip = f and t = t = 

(the dotted lines). From these results, we can see that, 
the superfluid phase of flat-band model becomes larger 
compared with the (/?-flux Bose-Hubbard model with only 
nearest hopping term. 



IV. COLLECTIVE MODES IN MOTT PHASE 



In this section, we first derive the effective action of 
this bosonic checkerboard model of flat-band. Then we 
calculate the dispersion and the energy gap of excitations 
in Mott phase. 



A. The effective action 

With the complex functions 6*(r) and bi(r) defined, 
the grand canonical partition function can be written as 



Z = Tre" 



36 = J Vb* 



vbe -S[b\b]/h 



(15) 



where Vb* and Vb denote functional integration for 6* 
and 6, respectively. The action S[b*,b] is given by 



S[b*,b] = / dr 
Jo 



d 



i ^ 
^£&*(t)&*(t)6;(t)Mt) 



(16) 



with /3 = &b the Boltzmann constant, and T the 

temperature. With a Hubbard- Stratonovich transforma- 
tion, we can rewrite the action as 



S[b*,b,r,i>] = S[b*,b] + / dr^(^(T)-6*(r))t lJ (^(r)-6 J (T)) 

= dT E 6 *( r )(^-^)^( r ) + ^E 6 *( r ) 6 *( r )^( r )^( T ) 

- J2 {K{r)U^j{r) + rdr)UMr)) + £ V>* (r)f^- (r) , 



(17) 



where -0* and -0 are the order parameter fields. Hnece where is the action with no hoping terms, 
we have Using the cumulate expansion formula 

e -s°«irM =exp|-i^rfrE^(r)^(r)J = e <*) + *«A? )-{A^) + ... (1Q) 

x y X>6*I>6exp ^ — 6] /ft — — dr we can e ^ ec ^ ve ac tion U P to second-order 

Xtii(6:(r)^-(r) +^?(r)6 J -(r)) I , (18) 

u / ' where perturbation term is 
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rnp / / pup 

Jo 



5(0) 



1 

~2h 



13 

h(3 r h(3 



^0 



drdr' Ujtvj. [b*(r)^ (r) + ^( T )b 3 (r)] [6*, (r')W (r') + V*< (r')&i' (r')] ) (21) 

iji'j' I 5(0) 

I 



With the correlations in the unperturbed system 

( b i b j) S W = ( b i b o)sW = 

(Wsw = (M*) 5(0) = ( b iK)sio) S id , (22) 

we get 

[</>*, iP}= dr^ ~T / drdr' 

Jo tj n Jo Jo 

E (T r [^(r)&*,(^)])5(o) • (23) 

iji'j' 

In general, the first term of above equation can be 
written as 

5>V*^ = £ ((^ k (r),^k(r))w(^[^ 

zj k 

(24) 

by a Fourier transformation, where T r is the imaginary- 
time order operator. And the quadratic term becomes 

E (lV[6r 1 (r)6; i ,(r')]) s(o) ^,(r') 



jii'j' 



E W (T T [6 ri (r)6 r * t (r')] ) s(0) V> v (r') 



= <T T [6 r , (r)&; ( (r')]) s(0) E (r)^., (r') . (25) 

It can be show that 

E^VV^)Vv.(t') 

ff 



I ^k(r') (26) 



Vsk(r') 



Hence at the transition point where 5^°) vanishes, the 
effective action can be written as 



This is consistent with the results of a two-site model in 
Ref.lnj. 



B. Dispersion of collective modes 

After obtaining the effective action of this bosonic 
checkerboard model of flat-band, we calculate the crit- 
ical point of MI-SF transition and then the dispersion of 
collective modes. 

For the bosonic checkerboard model of flat-band, the 
hoping terms are given by 



f. ■ — f*. 



( te i<p for NN terms, 
±t' for NNN terms, 
t" for NNNN terms, 

k otherwise. 



(28) 



Then the "H matrix was obtained asQ 

n 



A C* 
C B 



(T T [6(r)6*(r')])E(^k(r),^ k (r))^ 2 ' 



^sk(r') ; • 

(27) 



k k 

=4t" cos k x cos k y I + 4t cos (^(cos -y- cos ~^) <J x 
k k 

+ 4t sin <^(sin sin -^) cr y + 2^ (cos k x — cos k y )a z , 

(29) 

where / and cr x , a y , <j z is the identity and Pauli matrices. 
Hence 

A = 2t r (cos k x — cos k y ) + 4t" cos ^ cos fc^ , (30) 
23 = — 2t / (cos k x — cos k y ) + 4£" cos fc^ cos fc^ , (31) 

C = 4t cos cp cos — cos — + 4it sin <z? sin — sin — . 
^22 ^22 

(32) 

The (T T [b(T)b* (t 1 )]) S ( ) is just the negative of Matsub- 
ara Green's function and it can be evaluated as 

(T T [b(r)b*(T')}) s(0) = 6{t - r')(n + i) e -(-^)(— ')/« 

+ fl ( T '_ T ) ne (».-(n-l)tf)(r'-T)A _ 

(33) 
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After Fourier transformation to Hubbard- Stratonovich 
fields in the Matsubara frequencies, 



Bk(r) = -^J2 e ~ l0JT ^ B]i ^ ' 



(34) 



the effective action is given by 

s eS ir,4>] 

= E^km^Bk,™) (n - n 2 u m ) (f Bk >< 



k , rn 



k , rn 



(35) 



where 
and 



n + 1 



-ihuj m — fi + n/7 zfkj m + /i — (n — 1)C7 



fiG 1 (k,zo; m ) 

(.4 + B)C^ m -C (S 2 +C*C)/^-S 



(37) 



By substituting io; m — > cj m , we get the equation for real 
energy det[G _1 ] = 0, namely 

(AB - C*C) [(AB - C*C)fi m -(A + B)f Um + l] = 0. 

(38) 




FIG. 3: The energy gaps of a pair of quasi-particle and quasi- 
hole, AEi and AE2 via U when n = 1 and /x = 1. One can 
see that MFSF transition occurs at U = 12.15*. At MI-SF 
transition, Ai£i is still finite while AE2 close at the critical 
point. Hence the MI-SF transition is determined by AE2. 



Solve Eq ([36]) and Eq (|38j) . we can get the quasi- 
particle, quasi-hole spectra, 



Hui,qh 



fiW2,qh = \ 

where 



-2/x + (2n- 1)ET 
-2/x + (2n- 1)Z7 ■ 
-2/x + (2n- 1)C/ ■ 

-2/x + (2n- 1)ET 



1 



f 2 



f 2 



f 2 

/ com- 



f 2 

J com- 



A£7i >k 
AEi, k 
A£ 2 ,k 

A£ 2 , k 



(39) 



/[/ 2 



(4n + 2)C7 



/c 



u;ra+ 



f 2 

(40) 



, (36) A£; 2 ,k = fiw 2)9 p 



IIP 



(4n + 2)U 

fcurn — 



.4 + g ± y/(g - ff) 2 + 41^ 

2(„4£-C*C) 



(41) 
(42) 



Hence AEi/ 2 ,k are the dispersion of elementary excita- 
tions - a pair of quasi-particle and quasi-hole. 





FIG. 4: The spectra of quasi-particle and quasi-hole for the 
case of critical point, U « 12.15 \t\ when n = 1 and \i = 
1. The upper one is for huoi,qp, q h while the lower one is for 

fyjJ2,qp,qh- 

From the information of quasi-particle-quasi-hole spec- 
tra, we can determine the phase boundary of MI-SF tran- 
sition. We plot the energy gaps via U in Fi^3]with suit- 
able k where the energy gaps take their maximum value, 
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FIG. 5: The dispersion of a pair of quasi-particle and quasi- 
hole A£?2,k at the critical point with U ~ 12.15 \t\ when n — 1 
and fi=l. 



AE X 



/2 



The MI-SF transition oc- 



curs at/7 = ?7 c wl2.15|t| when the energy gap vanishes, 
which implies that A = AE2 = 0. From the results, 
we found that the topologically nontrivial flat-band only 
lightly changes the phase boundary between superfluid 
phase and Mott phase. 

In FigHJ we plot the spectra of two branches of quasi- 
particle and quasi- hole, hwi^ qp ^ q h and huj2,q P ,qh- One can 
see that in MI phase, the dispersion of a pair of quasi- 
particle and quasi- hole AEi^ = huj\ AV — fajj\, q h are al- 
ways have larger energy than Ai^,k = f^2,q P — ^2,qh- 
Thus we focus on the quasi-particle-quasi-hole excita- 
tions with dispersion A£?2,k in the following parts. For 
this case (U = U c ~ 12.15 the dispersion relation of 
the two quasi-particle quasi- hole spectra A£?2,k is dis- 
played in Figl5l At MI-SF transition we found that 
the dispersion of quasi-particles and quasi-holes A£?2,k 
show nodal-like behavior near special points in momen- 
tum space at k = (0, ±7r/2), k = (±7r/2, 0). 

In MI region, near MI-SF transition we also plot the 
spectra of two branches of quasi-particle and quasi-hole, 
f^i,qp,qh and hoj2,qp,qh fc> r the case of U = 12.5 \t\ in Figj6j 
In addition, the dispersion of a pair of quasi-particle and 
quasi-hole is shown in FigJTJ from which we found that 
the quasi-particle-quasi-hole excitations have energy gap 
and the dispersion of a pair of quasi-particle and quasi- 
hole A£?2,k becomes flat. To illustrate this effect, we 
calculate the flatness ratio, p = A/W where W is the 
band width of the energy spectra and A is the energy gap 
of of the excitations. The flatness ratio p for U = 12.5 \t\ 
is about 1.8. 

In MI region, far from MI-SF transition, we plot the 
spectra of two branches of quasi-particle and quasi-hole 
f^i,qp,qh an d J^2,qp,qh and the dispersion for this case 
of U = 16 \t\ in FiglH] and Figj9j respectively. In this 
figure, we can see that there exists a flat-band of A£?2,k 




FIG. 6: The spectra of quasi-particle and quasi-hole for the 
case of U = 12.5 \t\ when n = 1 and ji = 1. The upper one is 
for huJi, qPi qh while the lower one is for huo2,qp, q h- 



obviously. Now the flatness ratio p for U = 12.5 \t\ is 
about 16. 
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FIG. 7: The dispersion of the quasi-particle quasi-hole A^ 2 ,k 
for U — 12.5 \ t\ when n = 1 and ji = 1. 



Above results indicate that the dispersion of a pair 
of quasi-particle and quasi-hole will become more and 
more flat when we increase interaction strength, U. In 
Fi jlQ| we displayed the inverse of flatness ratio p via U 
of generalized Bose-Hubbard on a checkerboard model 
with topologically nontrivial flat-band. From these re- 
sults, we can see that the flatness ratio of the generalized 
Bose-Hubbard on a checkerboard model with topologi- 
cally nontrivial flat-band increases with increasing of the 
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FIG. 8: The spectra of quasi-particle and quasi- hole for the 
case of U = 16 \t\ when n = 1 and fi = 1. The upper one is 
for ftwi, qv , q h while the lower one is for ftw2, qP , q h- 




3 2 10 12 3 

FIG. 9: The dispersion of a pair of quasi-particle and quasi- 
hole Ai?2,k for U = 16 \t\ when n = 1 and ji = 1. 

interaction strength. On the other hand, we also calcu- 
late the flatness ratio of traditional Bose-Hubbard model. 
See blue line in FiglTOl From it we can see that the 
flatness ratio of traditional Bose-Hubbard model changes 
much more slowly than the flat-band model with increas- 
ing of the interaction strength. In this sense we can 
say that in MI phase, there indeed exist flat bands for 
(bosonic) quasi-particle or quasi-hole for the generalized 
Bose-Hubbard on a checkerboard model with topologi- 



cally nontrivial flat-band. 




FIG. 10: The inverse of flatness ratio 1/p via U / \ t\ when n = 
1 and fi = 1, where the red line is that of flat-band model while 
blue line is that of traditional Bose-Hubbard model. One can 
see that flatness ratio of the generalized Bose-Hubbard on 
a checkerboard model with topologically nontrivial flat-band 
increases much more rapidly, implying that there exists Hat- 
band for the bosonic excitations. 



V. CONCLUSION 

In this paper, using a decoupling approximation, we 
studied the generalized Bose-Hubbard on a checker- 
board model with topologically nontrivial flat-band in 
the mean-field level. We find that the MI-SF phase tran- 
sition of the flat-band case only lightly changes compared 
with the traditional Bose-Hubbard model. We also cal- 
culate dispersion relations of the collective modes Mott 
phase. The results show that in MI phase the (bosonic) 
quasi-particle or quasi-hole also has flat bands. In the 
end, we should point out that until now we have no idea 
about how to realize this bosonic checkerboard model 
of flat-band in optical lattice of cold atoms. In the fu- 
ture we will revisit this issue and find way to realize this 
bosonic checkerboard model of flat-band with particular 
big nearest-neighbor (NN) and the next-nearest-neighbor 
(NNN) hoppings. 
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